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ABSTRACT 

We present the results of an effort to identify and study a sample of the likely pro- 
genitors of elliptical (E) and lenticular (SO) K+A galaxies. To achieve this, we have 
searched a sample ^11,000 nearby (m r <16) early-type galaxies selected by morphol- 
ogy from the Sloan Digital Sky Survey (SDSS) Main spectroscopic sample for actively 
star-forming E and SO galaxies. Using emission line ratios and visual inspection of 
SDSS g-band images, we have identified 335 galaxies from the SDSS Fourth Data Re- 
lease (DR4) as actively star-forming E and SO galaxies. These galaxies make up about 
3% of the total early-type sample and less than 1% of all Main galaxies with m r <16. 
We also identified a sample of ^400 K+A galaxies from DR4 with m r <16; more than 
half of these are E and SO galaxies. We find that star-forming early-type galaxies and 
K+A galaxies have similar mass distributions; they are on average less massive than 
typical early-type galaxies but more massive than the average star-forming galaxy. 
Both of these types of galaxies are found in higher fractions among all galaxies in 
lower density environments. The fractions of star-forming E and SO galaxies and E 
and SO K+A galaxies depend on environment in nearly the same way. Model spectra 
fit to the stellar continua of the star-forming E and SO galaxies showed that their 
properties are consistent with star formation episodes of <1 Gyr in duration. The 
modelling results imply that on average, the star formation episodes will increase the 
stellar masses by about 4%. The results also imply that only episodes that increase 
the stellar mass by more than 2-5% will lead to K+A galaxies as we have defined them 
and that this is true for roughly 80% of the star forming E and SO galaxies in our 
sample. The estimated typical increase in stellar mass implies that new stellar com- 
ponents of about 2 x 10 8 M Q will be formed on average. There is also evidence that 
the star-forming regions within these galaxies are rotationally supported. These two 
results, when combined with the galaxies' total masses and lack of prominent disks, 
suggest that the star formation within these galaxies may be confined to relatively 
small, central disks, similar to the nuclear stellar and dust disks found in many low 
mass early-type galaxies. 

Key words: galaxies: elliptical and lenticular, cD - galaxies: star-burst - galaxies: 
stellar content 



1 INTRODUCTION 

The formation and evolution of early-type (elliptical and 
lenticular) galaxies is crucial to understanding galaxy forma- 
tion. The fact that the stellar populations in typical early- 
type galaxies seem to have formed long ago (~10 Gyr; e.g., 



* E-mail: helmbold@unm.edu 



Bower et al. 1992; Terlevich et al. 2001; Thomas et al. 
2005) appears to support the hypothesis that these galax- 
ies formed relatively quickly at high redshift and have had 
little or no new star formation since then. However, there 
is evidence that the stellar populations of early-type field 
galaxies may have formed at lower redshifts (van Dokkum & 
Franx 2001; Treu et al. 2002). There are even some early- 
type galaxies that show evidence of star formation within 
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the last ~1 Gyr (Zabludoff et al. 1996; Goto et al. 2003; 
Quintero et al. 2004; Yi et al. 2005; Kaviraj et al. 2007). 
These galaxies belong to a larger class of galaxies referred to 
as K+A (or E+A 1 ) galaxies discovered by Dressier & Gunn 
(1983). K+A galaxies are defined by their relatively strong 
stellar Balmer absorption lines and lack of nebular emission 
lines (usually Ha and/or [Oil]). These characteristics imply 
that while these galaxies formed some of their stars within 
the last ~1 Gyr, they are not currently forming stars. Recent 
studies of relatively large samples of K+A galaxies (Quin- 
tero et al. 2004; Goto et al. 2003) have revealed that these 
galaxies are preferentially found outside of galaxy clusters, 
implying that they may provide useful information regarding 
the evolution of early-type field galaxies. A detailed study of 
the environments of 266 K+A galaxies performed by Goto 
(2005) has shown that the local galaxy density within 100 
kpc of these galaxies is larger than what is found for other 
galaxies; no such excess was found on size scales typical of 
galaxy clusters or large scale structure. This implies that 
these galaxies may also be tracers of galaxy-galaxy tidal in- 
teractions or possibly merger/accretion events in the field. 

Understanding the role K+A galaxies play within the 
domain of galaxy evolution depends highly on our knowl- 
edge of the mechanisms which trigger and halt the star for- 
mation episodes that lead to these objects. To constrain the 
possible variety of such mechanisms and to obtain a more 
complete picture of these objects, a sample of likely progen- 
itors of K+A galaxies is required. While most K+A galax- 
ies are early-type galaxies, many also have prominent disk 
components and spiral structure (Goto 2005). Identifying 
the progenitors of these K+A spiral galaxies may be dif- 
ficult since current star formation is commonplace among 
spiral galaxies. A K+A spiral may be the result of the ces- 
sation of star formation within a relatively normal spiral 
as its gas is stripped via tidal interactions or ram pressure 
stripping within a group or cluster of galaxies (e.g., Verdes- 
Montenegro et al. 2001; Vogt et al. 2004). A K+A spiral 
may also be produced by a relatively short episode of star 
formation contained within the nucleus of a bulge domi- 
nated spiral. Conversely, current star formation is quite rare 
among early-type galaxies, and it is extremely unlikely that 
a process such as gas stripping could produce an early-type 
K+A galaxy. While some have argued that mergers of spiral 
galaxies could produce early-type K+A galaxies (Bekki et 
al. 2005), the optical and near-IR properties of the major- 
ity of K+A galaxies are consistent with the K+A spectral 
signature of these galaxies resulting from a short episode of 
star formation superimposed on an older stellar population 
(Balogh et al. 2005). It therefore seems promising to try to 
cull a sample of candidate K+A progenitors by identifying 
early-type galaxies that are actively forming stars. 

Recently, Fukugita et al. (2004) have discovered two 
galaxies in the Sloan Digital Sky Survey (SDSS) that arc 
unambiguously elliptical galaxies that are currently form- 
ing stars. To investigate the possibility that star-forming 
early-type galaxies like the elliptical galaxies discovered by 
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Fukugita et al. (2004) are plausibly the progenitors of early- 
type K+A galaxies, we have selected a large sample of ac- 
tively star-forming early-type galaxies using the SDSS spec- 
troscopic data and images. The sample selection, including 
the definition of different galaxy categories, is detailed in §2. 
The general properties of the sample galaxies and how they 
compare to K+A galaxies and early-type galaxies in general 
is discussed in §3. The distribution of star formation time 
scales among star forming early-type galaxies is explored us- 
ing population synthesis models. The model implementation 
and results are discussed in §4; further applications of the 
model results are discussed in §5. In §6, the main results are 
summarised and useful follow-up observations of this sample 
are discussed. 



2 SAMPLE SELECTION AND GALAXY 
CATEGORIES 

2.1 Identifying star-forming galaxies 

The selection of our sample of actively star-forming early- 
type galaxies consists of two basic parts. The first involves 
the selection of a large sample of actively star-forming galax- 
ies. The second consists of the combined use of morphologi- 
cal selection criteria and visual inspection of optical images 
to produce from this sample a sub-sample dominated by el- 
liptical and lenticular galaxies. We stress that the selection 
process was designed to yield a relatively clean sample of 
star forming elliptical and lenticular galaxies and is some- 
what complex. This makes a meaningful comparison of the 
number of galaxies in the resulting sample to those in other 
samples difficult, and such comparisons are not part of the 
analysis and interpretation that follows. The two parts of 
this selection process are detailed below. 

All of the galaxy data used here are taken from the 
fourth data release (DR4) of the SDSS (Adelman-McCarthy 
et al. 2006). The SDSS is an imaging an spectroscopic sur- 
vey of 8,000 square degrees of the northern sky. From ugriz 
images, ~ 10 6 spectroscopic targets were selected and ob- 
served with a fibre fed spectrograph that can observe 640 
objects simultaneously (see Gunn et al. 1998; York et al. 
2000; Stoughton et al. 2002). These targets were selected 
in an automated way (Eisenstein et al. 2001; Strauss et al. 
2002; Richards et al. 2002) based on photometric proper- 
ties measured within an automated data processing pipeline 
that performs astrometry (Pier et al. 2003) , source identifi- 
cation, de-blending, photometry (Lupton et al. 2001), and 
calibration (Fukugita et al. 1996; Hogg et al. 2001; Smith 
et al. 2002; Ivezic et al. 2004) . The fibre placement was also 
automatically determined for the targeted objects (Blanton 
et al. 2003a). 

We have identified star-forming galaxies within the flux 
limited Main sample (Strauss et al. 2002) from DR4. While 
the Main sample has a flux limit of m r < 17.77, we have 
used a limiting magnitude of m r < 16 for reasons explained 
below. To identify actively star forming galaxies, we have 
used emission line ratios and a diagnostic diagram similar 
to those used by Osterbrock & Ferland (2006). For this pur- 
pose, we have used the emission line fluxes measured from 
SDSS spectra by Tremonti et al. (2004). The basic proce- 
dure used by Tremonti et al. (2004) is as follows. The fluxes 
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were measured by first fitting to each galaxy spectrum a lin- 
ear combination of model spectra with different star forma- 
tion histories produced using the population synthesis code 
of Bruzual & Chariot (2003). This fit was then subtracted 
to remove the stellar continuum, a step which is vital to our 
efforts to obtain accurate emission line fluxes for early-type 
galaxies whose stellar continua, especially near the galax- 
ies' centres, are typically quite bright. After subtraction of 
the model continuum fits, Gaussian functions were fit to the 
galaxy's emission lines. From these fits, the flux, equivalent 
width, and velocity width was obtained for each emission 
line as well as estimates for the errors in each of these quan- 
tities. 

To complete the first step in our selection process, we 
have chosen to use the emission line ratios [OIII]/H/3 and 
[Nil] /Ha to select actively star-forming galaxies and exclude 
Seyfert and LINER galaxies. We have chosen to only include 
galaxies with m r <I6 that are nearby enough that any spi- 
ral arms or prominent disks will be apparent in their g-band 
images. We have chosen this flux limit because it similar 
to the limit used by Nakamura et al. (2003) for their cata- 
logue of morphologically classified SDSS galaxies. They have 
demonstrated that from the SDSS images, it is possible to 
obtain reasonably accurate morphological classifications by 
eye for galaxies down to this magnitude limit. Formally, the 
limiting r'-band magnitude used by Nakamura et al. (2003) 
was 15.9. However, the fluxes they used were obtained us- 
ing an older version of the SDSS photometric pipeline using 
the photometric system of SDSS early data release (EDR; 
Stoughton et al. 2002), the u 1 , g 1 ,r' ,i' , z system. In the 
photometric system used in subsequent data releases (i.e., 
u, g,r,i, z), including DR4, 4-5% of the galaxies from their 
sample have m r >15.9. In light of this, we have chosen to 
increase the limiting magnitude to 16. 

In Fig. 1, we plot the two emission line ratios for all 
galaxies with m r <16 and > 5a detections (i.e., the emis- 
sion line fluxes are at least five times the uncertainty given 
by Tremonti et al. (2004)) of the [NII]A6583, H/3, and Ha 
emission lines. In the upper panel of the figure, we only in- 
clude those galaxies with > 3<r detections of the typically 
weaker [OIII]A5007 line. For galaxies that have no signifi- 
cant [OIII] emission but have > 5a detections of the other 
three emission lines, we plot upper limits for [OIIIJ/H/3 in 
the lower panel of Fig. 1. We also plot the theoretical upper 
limit for star-burst galaxies from Kewley et al. (2001) as 
a dashed curve in both panels of Fig. 1. From these plots, 
it is apparent that the star-forming and AGN galaxies form 
two relatively distinct sequences and that even though they 
are consistent with star-burst galaxies, many of the galaxies 
that lie below the Kewley et al. (2001) boundary may con- 
tain significant AGN components. Since the vast majority 
of elliptical and lenticular galaxies have emission line ratios 
consistent with extended LINER-like emission (Phillips et 
al. 1986; Goudfrooij et al. 1994; Sarzi et al. 2006), we do 
not want to use such a generous upper limit for the selection 
of our sample since our ultimate goal is to identify elliptical 
and lenticular galaxies that are most likely forming stars. 
We have therefore adjusted the curve from Kewley et al. 
(2001) to provided a more conservative upper limit which is 
given by 



( [OIII]A5007\ 



0.61 



log ([Nil] /Ha) -0.12 



+ 1.14 



(1) 



This upper limit is plotted as a solid curve in both panels of 
Fig. 1 and is similar to other empirically derived limits such 
as that used by Kauffmann et al. (2003a, see Fig. 1). 



2.2 Identifying early-type galaxies 

With the boundary given by equation (1), we can select a 
sample of star-forming galaxies within which we can search 
for elliptical and lenticular galaxies. Using this boundary, we 
find ~11,000 star-forming galaxies with m r <16. Morpho- 
logically classifying such a large sample by eye is impracti- 
cal. Therefore, we must isolate those star-forming galaxies 
that are most likely to be early-type galaxies so that we 
may obtain a more manageable number of galaxies for vi- 
sual inspection to eliminate spiral contamination. We have 
chosen to do this using a combination of the Sersic indexes, 
n s , and the axis ratios, b/a, for the galaxies. The Sersic in- 
dexes are taken from the Sersic fits performed by Blanton et 
al. (2003b) to surface brightness profiles measured from the 
i-band images by the SDSS photometric pipeline using cir- 
cular apertures. Tests performed by Blanton et al. (2003b) 
indicate that these fits perform reasonably well down to 
the magnitude limit of the Main sample for galaxies with 
6/o^0.5. For the quantity b/a, we used values from the best 
fitting de Vaucouleurs profile as determined by the SDSS 
photometric pipeline using the i-band images which take 
the size of the seeing disk into account (see Stoughton et al. 
2002). The fits are performed assuming that the isophotal 
axis ratio does not change with radius, and consequently, 
one value of b/a is solved for for each galaxy. To determine 
the best combination of these two parameters, we have used 
the catalogue of 1875 morphologically classified SDSS galax- 
ies compiled by Nakamura et al. (2003). We adjusted the 
imposed limits for n s and b/a so that the largest number of 
the elliptical galaxies in the Nakamura ct al. (2003) cata- 
logue would be retained while at the same time, the largest 
number of galaxies of type Sa or later would be rejected. 
The criteria we derived were n s >3.5 and b/a >0.7. From 
the Nakamura et al. (2003) catalogue, 67% of the galaxies 
that pass these two cuts are E or E/S0 galaxies, 16% are SO 
or SO/a galaxies, 16% are of type Sa or later, and 1% are 
irregular or peculiar galaxies. About 90% of all the elliptical 
galaxies within the Nakamura et al. (2003) catalogue meet 
the two criteria. We note that 25% of the galaxies within the 
Nakamura et al. (2003) catalogue that have n a >3.5 are spi- 
ral galaxies, implying that the additional cut based on axis 
ratio is useful. We also note that since the Sersic fits are more 
reliable for galaxies with larger axis ratios and that the b/a 
values taken from the SDSS pipeline are more appropriate 
for galaxies with n s ~4, the two criteria we have chosen 
imply that the measured values for n s and b/a are reliable. 
Hereafter, when we refer to galaxies as "early- type", we are 
referring to those galaxies with n s >3.5 and b/a >0.7; those 
star-forming galaxies that meet these criteria are referred to 
as star- forming early-type, or "SFE" galaxies. 

Among the star-forming galaxies from DR4 that have 
m r <16, 752 have n 3 >3.5 and b/a >0.7, a reasonable num- 
ber of SFE galaxies for which visual inspection of the images 
can be done in a practical amount of time. We have exam- 
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Figure 1. For galaxies with m r <16 and > 5<r detections of the [NIIJA6583, H/3, and Ha emission lines, [OIII]/H/3 versus [Nil] /Ha for 
galaxies with > 3c detections of the [OIII]A5007 line (upper panel) and upper limits for [OIIIJ/H/3 versus [Nil] /Ha for galaxies with 
no significant detection of [OIII] emission (lower panel). In both panels, the dashed curve is the upper limit for star-burst galaxies from 
Kewley et al. (2001), the dotted line is the empirical upper limit for star forming galaxies used by Kauffmann et al. (2003a), and the 
solid curve is the more conservative limit used here. Early-type galaxies (see §2) are represented by X 's; star forming early-type (SFE) 
galaxies are represented by white circles for elliptical and lenticular (E) galaxies and white triangles for early-type spiral (S) galaxies. 



ined the p-band images of these galaxies and have put them 
into three categories, (i) elliptical and lenticular galaxies, (ii) 
galaxies that have prominent disk components, spiral arms, 
or disks with "lumpy" structure indicative of Hn regions 
(i.e., likely spiral galaxies), and (hi) peculiar, irregular, or 
merging galaxies as well as Hn regions within nearby galax- 
ies incorrectly identified as individual galaxies by the SDSS 
pipeline. For convenience, we hereafter refer to these three 
groups as E, S, and P. The numbers of SFE galaxies within 
these groups are summarised in Table 1. The p-band images 
and spectra of four SFE E galaxies and four SFE S galaxies 
chosen at four different r-band magnitudes (m r ~12.5, 13.5, 
14.5, and 15.5) are displayed in Fig. 2. It should be noted 
that we have not performed a detailed morphological classi- 



fication of our SFE galaxies. For instance, all galaxies with 
dominant spheroidal components and no prominent disks or 
spiral structures were placed within the E category, regard- 
less of any slightly atypical features such as off-centre nuclei 
(e.g., see the g-band image of SDSS J121758.08+071626.22 
in Fig. 2), and galaxies within the S category were not bro- 
ken up into subclasses such as Sa, Sb, Sc, etc. 

While Fig. 2 demonstrates that our visual inspection 
has worked well for those cases, it does little to reassure us 
that our morphological classification is consistent over the 
entire range in redshifts / apparent size inhabited by our sam- 
ple galaxies. We can roughly test the quality /reliability of 
our visual classifications by comparing our visual inspection 
as a function of apparent magnitude to that of Nakamura 
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Figure 2. Spectra (log versus log A) from 3800 to 8000 A and 3-band images of four SFE E galaxies (left, see §2) and four SFE S 
galaxies (right) chosen at four different r-band magnitudes (m r ~12.5, 13.5, 14.5, and 15.5). 



et al. (2003) to at least see if our classifications are as con- 
sistent as theirs are as a function of distance. To do this, 
we have computed the cumulative distributions of m r for 
all elliptical and lenticular galaxies (type SO/a or earlier) 
from the Nakamura et al. (2003) sample that have n s >3.5 
and b/a >0.7 and for all spiral galaxies (type Sa or later) 
that also have n s >3.5 and b/a >0.7. These distributions 
are plotted in Fig. 3 along with the same distributions for 
our SFE E and SFE S galaxies. Each distribution has been 
normalised by the number of galaxies with m r <15.9 since 
the Nakamura et al. (2003) sample is significantly incom- 



plete for magnitudes fainter than this limit. From Fig. 3, it is 
apparent that the shapes of the distributions for our SFE E 
and SFE S galaxies match those of the corresponding Naka- 
mura et al. (2003) sub-samples. This implies that while we 
have no direct test of how well our morphological classifica- 
tion works with redshift, we can at least be sure that our 
visual inspection has performed as well as that of Nakamura 
et al. (2003) as a function of apparent magnitude. 

Since we are reasonably confident in our morphological 
classifications, we have used these classifications to eliminate 
spiral galaxies (i.e., SFE S galaxies) from our analysis. As 
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Figure 3. A comparison of the visual classification of our SFE 
galaxies to that of the Nakamura ct al. (2003) sample of SDSS 
galaxies. The cumulative distributions of the r-band magnitude 
for E, E/SO, SO, and SO/a galaxies from the Nakamura ct al. 
(2003) catalogue that have n s >3.5 and b/a >0.7 is represented 
by the lower solid line. The distribution for the Sa or later-type 
galaxies from the Nakamura ct al. (2003) catalogue that have 
n s >3.5 and b/a >0.7 is plotted as the upper solid line. Simi- 
lar cumulative distributions are also plotted for SFE E (see §2) 
galaxies (lower dashed line), and for SFE S galaxies (upper dashed 
line). Each curve has been normalised to the number of galaxies 
at m r =15.9, which is similar to the r'-band (i.e., in the older 
SDSS photometric system) magnitude limit of the Nakamura et 
al. (2003) catalogue. The fact that the shapes of the curves for 
our SFE E and SFE S samples match those of the curves for the 
corresponding Nakamura et al. (2003) sub-samples implies that 
our visual inspection has performed as well as that of Nakamura 
et al. (2003) as a function of apparent magnitude. 



discussed in §1, it is more likely that elliptical and lentic- 
ular K+A galaxies result from episodes of star formation 
within morphologically similar galaxies while K+A spiral 
galaxies may result from either a similar scenario or from a 
truncation of star formation by processes such as gas strip- 
ping. This implies that while the SFE E galaxies are excel- 
lent candidates for the progenitors of morphologically sim- 
ilar K+A galaxies, it is substantially more difficult to tie 
the SFE S galaxies to spiral K+A galaxies in a similar way. 
This is made even more difficult by the fact that the 3 arcsec 
SDSS spectrograph aperture may miss significant amounts 
of disk star formation in nearby galaxies (see Fig. 2), imply- 
ing that many spiral galaxies may be misclassified as K+A 
using their SDSS spectra. For these reasons, the remainder 
of this paper will focus on the SFE E galaxies and how they 
relate to morphologically similar K+A galaxies. 

To roughly compare the range of star formation rates 
(SFRs) occupied by our SFE E galaxies with that for "nor- 
mal" star forming galaxies, we have plotted the Ha emission 
line equivalent width distributions for these galaxies in Fig. 
4. Since the Ha luminosity is a direct measure of the SFR 



Table 1. Galaxy Categories 



Category Number 



All SDSS galaxies with m r <16 39,621 

Early-type with m r <16 10,870 

S/N>5 for H/3, Ha, and [Nil] 16,282 
S/N>5 for H/3, Ha, [Nil], and S/N>3 for [OIII] 15,091 

AGN with m r <16 1,754 

Star Forming with m r <16 11,417 

SFE 752 

SFE E 335 

SFE S 402 

SFE P 15 

K+A, m r <16 435 

K+A, m r <16, S/N>3 for Ha 220 

K+A E 253 

K+A S 142 

K+A P 40 



(Kennicutt 1998), the Ha equivalent width, W(Ha), is a 
good tracer of the SFR per unit stellar mass. In Fig. 4, 
we have plotted separately the W(Ha) distributions for all 
star forming galaxies with m r < 16 and for SFE E galaxies. 
From these distributions, we find that in terms of W(Ha), 
our SFE E galaxies are relatively similar to more typical star 
forming galaxies with slightly (a factor of ~1.1 times) lower 
values of W(Ha) on average. For comparison, we have also 
plotted the W(Ha) distribution for AGN galaxies, which 
we take to be all galaxies with m r <16 that have values of 
[OIIIJ/H/3 that are greater than the Kewley et al. (2001) 
upper limit for star-burst galaxies (see Table 1) displayed 
in Fig. 1. The typical value of W(Ha) for these galaxies 
is factor of ~3 lower than what was found for star form- 
ing galaxies, significantly lower than what was found for 
SFE E galaxies. This implies that in terms of their levels of 
Ha emission, SFE E galaxies are more consistent with typ- 
ical star forming galaxies than with AGN. It should also be 
noted that the Ha equivalent widths for the SFE E galaxies 
are substantially higher than what can be produced by pho- 
toionisation by post-asymptotic giant branch stars (PAGB) 
within post-star-burst or older, quiescent stellar populations 
which ranges from about 0.3 to 3 A according to the mod- 
els of Binette et al. (1994). This is in contrast with typical 
early-type galaxies where the level of Ha emission is usually 
consistent with both LINER galaxies and model predictions 
of photoionisation by PAGB stars (e.g., Goudfrooij et al. 
1994; Binette et al. 1994; Macchetto et al. 1996). 

2.3 Identifying K+A galaxies 

Since our major goal is to explore the possibility that our 
SFE E galaxies will become morphologically similar K+A 
galaxies, we need to have a sample of these objects as well. 
Among all Main galaxies with m r < 16, we have identified 
K+A galaxies of all morphological types and have used vi- 
sual inspection of the g-band images to identify E and SO 
K+A galaxies. The first step in this process was to choose a 
definition for K+A galaxies. We have adopted a K+A def- 
inition that is very similar to that used by Quintero et al. 
(2004) to identify K+A galaxies using SDSS spectra. The 
main difference is that we have chosen to use the strength 
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of the H<5 absorption line instead of the quantity A/K used 
by Quintero et al. (2004). The quantity A/K is essentially 
the ratio of the flux between 3800 and 5400 A due to the 
recently (within ~1 Gyr) formed stellar component (i.e., the 
component that resembles a main sequence A star) to that 
due to the more quiescent component (i.e., the component 
that resembles a K giant star). We have opted to use the 
strength of the H<5 absorption line instead of A/K because 
the H8 absorption equivalent width is relatively insensitive 
to dust extinction when compared to the ratio A/K which 
is based on the shape of a galaxy's optical spectrum over a 
range of 1600 A in wavelength. 

To develop our own similar K+A definition, we have 
chosen to use the spectral index H<5a (Worthey & Ottaviani 
1997) and the Ha emission line equivalent width, W(Ha). 
The index H5a is designed to provide a measure of the 
equivalent width of the stellar H<5 absorption line and is 
defined such that positive values imply absorption. For the 
nebular Ha emission line equivalent width, W(Ha), we have 
chosen to associate positive values with emission because (i) 
such equivalent widths typically occupy a large range in val- 
ues and it is more convenient to work with values of log 
W(Ha) and (ii) we wish to be able to compare our efforts 
with similar work performed by previous authors who have 
used the same convention (e.g., Quintero et al. 2004). For 
both H8a and W(Ha), we used the values measured by 
Tremonti et al. (2004) where the values of H<5a were cor- 
rected by Tremonti et al. (2004) for emission lines that 
originate from the individual galaxies and in some cases, 
from foreground sky emission. In the upper panel of Fig. 5, 
we plot W(Ha) versus H<5a for galaxies with m r <16 and 
> 3a detections of the Ha emission line. For galaxies with 
no significant Ha detection, we plot upper limits for W(Ha) 
versus ~HSa in the lower panel of Fig. 5. In developing a 
definition for K+A galaxies, we sought to exclude actively 
star-forming and quiescent galaxies while including as many 
K+A galaxies as possible. To exclude star-forming galaxies, 
we have empirically derived an upper limit for W(Ha) as a 
function of Hoa using the lowest observed values for W(Ha) 
for star- forming galaxies as they are defined in Fig. 1. This is 
very similar to what was done by Quintero et al. (2004) who 
used an upper limit for W(Ha) which varied linearly with 
A/K. Our adopted upper limit is plotted as a solid line in 
both panels of Fig. 5. Formally, our K+A definition requires 
log W{Ha) < O.IIHSa + 0.15 and H<5 A > 2 A; this second 
requirement is designed to exclude quiescent galaxies. 

Using the K+A definition illustrated in Fig. 5, we found 
435 K+A galaxies within DR4 with m r <16 which make up 
~1% of all DR4 galaxies with m r <16. This is very similar to 
the K+A galaxies within the Quintero et al. (2004) sample 
which make up roughly 1% of galaxies from the SDSS Main 
sample. We note that the relative number of galaxies in ei- 
ther our K+A sample or the Quintero et al. (2004) sample is 
approximately an order of magnitude higher than has been 
found using more strict definitions of K+A galaxies (e.g., 
Zabludoff et al. 1996; Goto et al. 2003). We also note that 
roughly half of our K+A galaxies have > 3a detections of Ha 
(see Table 1 and Fig. 5) while only about 6% have > 3a de- 
tections of H<5 emission. We further note that both our K+A 
selection criteria and those of Quintero et al. (2004) do not 
include any requirement based on the flux of the [OII]A3727 
doublet, unlike some previous authors (e.g., Zabludoff et al. 
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Figure 4. The distributions for the Ho emission line equiva- 
lent width, W(Ha), for all star forming galaxies with m r <16 
(dashed line) and for SFE E galaxies (solid line). For comparison, 
the W(Ha) distribution for AGN galaxies with m r <16 (dotted 
line) is also plotted where we consider all galaxies with values 
of [OIIIJ/H/3 greater than the upper limit for star- burst galax- 
ies given by Kewley et al. (2001) illustrated in Fig. 1. For each 
histogram, Ntot refers to the total number of galaxies within the 
corresponding sub-sample. 



1996; Goto et al. 2003). Several authors have noted a sig- 
nificant fraction of galaxies with strong Balmer absorption 
and weak or no Ha emission do have significant levels of 
[Oil] emission (e.g., Liu & Kennicutt 1995; Zabludoff et al. 
1996; Goto et al. 2003; Quintero et al. 2004). About 18% 
of our K+A galaxies have > 3a detections of [Oil] emission 
with a typical equivalent width of ~ 10 A. Among those 
K+A galaxies with detected Ha and [Oil] emission, the ra- 
tio of [Oil] /Ha is consistent with what has been found for 
more typical early-type galaxies (Goudfrooij et al. 1994; Bi- 
nette et al. 1994). This implies that the [Oil] emission is 
not indicative of residual star formation and is more consis- 
tent with LINER-like emission (e.g., Goudfrooij et al. 1994; 
Yang et al. 2006), possibly associated with photoionisation 
by PAGB stars (Binette et al. 1994; Goudfrooij 1999). 

To identify elliptical and lenticular K+A galaxies we 
have used visual inspection of g-band images. Since there are 
only 435 K+A galaxies with m r < 16, as opposed to 11,417 
star forming galaxies, it was not necessary to apply the same 
morphological selection criteria to provide a smaller, more 
manageable list of images to inspect. We have therefore not 
required that the K+A galaxies have n s > 3.5 and b/a > 0.7. 
After inspecting the g-band images of the K+A galaxies, we 
put them into the same three morphological categories as 
were used for the SFE galaxies, E, S, and P. The numbers of 
K+A galaxies within each of these groups are summarised 
in Table 1. These morphological classifications provided us 
with a sample of 253 elliptical and lenticular K+A galaxies 
with which we can compare our SFE E galaxies. 



8 J. Helmboldt et al. 




Figure 5. The Ha emission line equivalent width, W(Hct) (emission is positive; sec §2.3), versus the strength of the H5 absorption line, 
H&a (absorption is positive; see §2.3), for galaxies with > 3cr detections of Ha and m r <16 (upper). For galaxies with no significant 
detection of Ha, the upper limit for W(Ha) is plotted in the lower panel. The point styles and are the same as those used in Fig. 1 
with elliptical and lenticular K+A galaxies represented by black circles and spiral K+A galaxies represented by black triangles. In both 
panels, the adopted definition of K+A galaxies is illustrated by the solid lines (see §2.3). 



3 EARLY- TYPE STAR FORMING AND 
EARLY- TYPE K+A GALAXIES 

If our SFE E galaxies will in fact become K+A E galax- 
ies, then they must have a few basic properties in common 
with the galaxies within the K+A E sample. To check this, 
we first used the 1/V ma x method to estimate the number 
density of SFE E and K+A E galaxies with z>0.05 and 
m r >13 within different bins of absolute z-band (i.e., 
K-corrected to a redshift of 0.1) magnitude, Mo.i z , stellar 
mass, M», and stellar line-of-sight (LOS) velocity dispersion, 
a v ; the results are plotted in Fig. 6. For these calculations, 
and for all other distance dependent calculations through- 
out this paper, we assumed (fi m , £Ia, h)=(0.3,0.7,0.7). The 
lower redshift limit of 0.05 was used to minimise the effects 
of aperture bias resulting from the use of a fixed angular 



aperture (i.e., the SDSS spectrograph 3 arcsec fibre aper- 
ture; Gomez et al. 2003). The lower magnitude limit of 
13 was implemented to exclude the most nearby galaxies. 
This was done because nearby galaxies with bright cores 
will not be targeted for spectroscopy by the SDSS pipeline 
to prevent crosstalk among fibres (Strauss et al. 2002), 
implying that the 1/Vmax method will provide an insuf- 
ficient completeness correction for these objects. The val- 
ues for Mo.i z and M, were taken from (Kauffmann et al. 
2003b) and include their K-corrections and extinction cor- 
rections. To obtain measurements of a v , we have used the 
IDL routine vdispfit written by D. Schlegel. This routine de- 
termines the best-fitting velocity dispersion and the la error 
in that dispersion by cross-correlating each spectrum with 
several template spectra that have been broadened by vari- 
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ous Gaussian velocity distributions while masking regions of 
the spectrum that may contain emission lines. For each of 
the two groups, the distributions have been normalised by 
the total estimated number density. This was done because 
the relatively complex selection process used to cull each 
of the sub-samples precludes any meaningful comparison of 
the actual number densities between the two groups, and we 
are therefore only interested in comparing the shapes of the 
sub-samples' distributions. 

We have also computed the number density of all star- 
forming (see §2.1) and all early-type (see §2.2) galaxies from 
our parent sample with z>0.05 and 13< m r <16 within bins 
of Mo.i z , M„, and a v . So that the shapes of these distribu- 
tions can be compared with those of the SFE and K+A 
galaxies, we have plotted the results in Fig. 6 with each 
distribution normalised by its total estimated number den- 
sity. From Fig. 6, it is apparent that for all three quantities, 
the distribution functions for SFE E and K+A E galaxies 
both lie somewhere in between those for normal star-forming 
galaxies and those for typical early-type galaxies. This im- 
plies that each of these types of galaxies is on average more 
massive than the typical star-forming galaxy, but less mas- 
sive than the average early-type galaxy. The shapes of the 
distributions for all three properties appear essentially iden- 
tical for both SFE E and K+A E galaxies. This implies that 
the galaxies from both groups have about the same stellar 
and dynamical masses on average. 

The results above imply that our SFE E galaxies have 
roughly the same masses as K+A E galaxies, but are they 
located in the same environments? To address this question, 
we have quantified the environment around each galaxy by 
measuring the projected distance, R4, to the fourth near- 
est Main galaxy with a radial velocity separation <1200 
km s _1 and with an absolute r-band magnitude brighter 
than -20. Given the Main sample limiting magnitude of 
17.77, this allows us to measure R4 for all galaxies out to a 
redshift of 0.079. We have again used the 1/V ma x method 
to compute the number density of all star-forming, early- 
type, SFE E, and K+A E galaxies with 13< m r <16 and 
0.05<z<0.079, this time within different R 4 bins. In Fig. 7, 
we plot the fraction of the total number density of Main 
galaxies with 13< m r <16 and 0.05<z<0.079 contained 
within each galaxy category for each of these bins. As may 
be expected, the fraction of star-forming galaxies drops sub- 
stantially at lower values of R4 whereas the fraction of early- 
type galaxies does the exact opposite. The fraction of both 
SFE E and K+A E galaxies is higher for the largest R4 
bin, implying that in terms of environment, these galax- 
ies have more in common with typical star-forming galaxies 
than most early-type galaxies. This is consistent with what 
was found for the K+A galaxy samples of Hogg et al. (2005) 
and Goto (2005) as well as the somewhat more evolved 
"E+F" galaxies from the sample of Nolan et al. (2007). It 
should also be noted that the shapes of the trends between 
the fractions of SFE E and K+A E galaxies and R4 are 
quite similar, implying that the two types of galaxies typi- 
cally reside within similar environments. Combined with the 
findings of the authors cited above, this indicates that it is 
plausible that our SFE E galaxies lie along a common evo- 
lutionary sequence with both K+A and E+F galaxies. 
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Figure 6. Galaxy number density as a function of Mo.i z (the ab- 
solute 2-band magnitude K-corrected to z=0.1, the median red- 
shift of the Main sample; upper), stellar mass, M, (middle), and 
velocity dispersion, cr v (lower), for SFE E (open circles) and K+A 
E (filled circles). We also plot the same for all star-forming galax- 
ies (dashed line; see §2.1) and all early-type galaxies (solid line; 
see §2.2) with m r <16 so that the shapes of the SFE and K+A 
distributions may be compared to the shapes of these distribu- 
tions. In all panels, the number density in each bin is divided by 
the total number density since our relatively complex selection 
criteria preclude any meaningful comparison of absolute number 
densities among the sub-samples. 



4 MODELLING THE STELLAR CONTINUA 
OF STAR-FORMING ELLIPTICAL AND 
LENTICULAR GALAXIES 

4.1 Basic assumptions 

The fact that our SFE E galaxies make up a relatively small 
fraction (~3%) of all early-type galaxies implies that the cur- 
rent star formation within these galaxies is likely short-lived. 
Given their morphological similarity to the majority of K+A 
galaxies which are roughly just as rare among all early-type 
galaxies, it seems likely that SFE E galaxies will evolve into 
early-type K+A galaxies. The fact that the mass distribu- 
tions and environments of our samples of SFE E and K+A E 
galaxies are similar is consistent with this conclusion. In this 
section we will explore the likely evolutionary link between 
our SFE E galaxies and early-type K+A galaxies in more 
detail. We will discuss our estimates of the magnitude and 
time-scale of the star formation in the early-type galaxies, 
and the subsequent spectral evolution of the SFE E galaxies. 
We will do this from model fits to their stellar spectra which 
utilise the population synthesis models of Bruzual & Char- 
lot (2003). We will also discuss how these fits were used to 
determine under which conditions the galaxies will appear 
as K+A galaxies with spectroscopic properties as defined by 
us earlier (see §2.3 and Fig. 5). 

To construct models that can be fit to the SFE E galaxy 
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Figure 7. Within bins of the distance to the fourth nearest Main 
galaxy within 1200 km s _1 with M r <-20, R4, the fraction of the 
total number density of Main galaxies attributed to star-forming 
(stars), all early-type (x's), SFE E (open circles), and K+A E 
(filled circles) galaxies. All galaxies used have 13< m r <16 and 
0.05<z<0.079 (see §3). 



spectra, the main assumption that we have made is that the 
SFE E galaxies consist of two basic stellar components, an 
older pre-existing population and a new population that is 
currently being formed. In fact, Balogh et al. (2005) con- 
cluded that the stellar H5 absorption equivalent widths and 
near-IR colours of many K+A galaxies are best explained 
by star formation histories where relatively new (within the 
last ~1 Gyr) stars were formed among existing older stel- 
lar populations. We note that no explicit assumption was 
made regarding the formation mechanism (e.g., monolithic 
collapse or major mergers) or formation epoch for the as- 
sumed underlying older stellar population; we only assumed 
that the gross spectral properties of this population can be 
well approximated with an instantaneous burst model at a 
relatively large age. We have also made the simplifying as- 
sumption that the star formation within the SFE E galaxies 
exponentially decreases with time. Since the stellar continua 
alone cannot be used to effectively constrain the e-folding 
time, t, for this scenario, we have chosen to perform four 
separate fits for each spectrum, assuming r = 100 Myr, 
300 Myr, 1 Gyr, and 00 (i.e., a constant SFR). To further 
simplify our modelling efforts, we have chosen to only use 
Bruzual & Chariot (2003) models that assume the initial 
mass function of Chabrier (2003). These assumptions allow 
us to use the Bruzual & Chariot (2003) instantaneous burst 
models to compute suites of model spectra for the older pop- 
ulation and the exponentially decreasing and constant SFR 
models to compute the same for the younger stellar com- 
ponents. To this end, we have constructed a suite of 1100 
spectra for each of these models using five different metallic- 
ities (Z=0.0004, 0.004, 0.008, 0.02, and 0.05) and 220 times 
steps ranging from 10 5 ' 1 to 10 10 ' 3 yr. Each of the model 



spectra was normalised by dividing by the total stellar mass 
at the appropriate time step which includes the Bruzual & 
Chariot (2003) computation for the amount of stellar mass 
loss. 



4.2 Parameter constraints 

Our basic approach for estimating the age of the younger 
population within each SFE E galaxy was to find the best 
fitting linear combination of an instantaneous past burst 
spectrum and one of the exponentially decreasing SFR spec- 
tra from the five suites of model spectra. The coefficients in 
this combination are essentially the stellar masses of the 
two component populations. However, in practice, we must 
recognise that (i) the stellar mass, age, and metallicity are 
not the only factors that affect the shapes of the spectra and 
(ii) the effects of different properties on the spectral shapes 
may be degenerate (e.g., age and metallicitiy) . The two ad- 
ditional properties that must be considered when finding the 
best fitting model spectrum are the amount of internal dust 
extinction and the velocity dispersion. To simplify matters, 
we assumed one velocity dispersion for the combined model 
spectrum, the dispersion obtained for the actual galaxy spec- 
trum using vdispfit (see §3). We have also made the simplify- 
ing assumption that the old and new populations suffer the 
same amount of dust extinction. This assumption is rela- 
tively good if the star forming regions occupy an area greater 
than or equal to the area covered by the SDSS spectrograph 
fibre aperture. Narrow-band imaging of a few or our SFE 
E galaxies suggests that this is indeed the case (Helmboldt 
2005). For each spectrum, we therefore must fit for several 
different parameters and we must constrain each of these as 
strictly as possible to overcome the degeneracy between the 
age of the younger population (the main parameter we are 
interested in) and any other parameter. We have done this 
in the following way: 

(i) Metallicity: We have good estimates for the oxygen 
abundance for the ionised gas within each SFE E galaxy 
from the work of Tremonti et al. (2004). It is likely that 
the ionised gas is associated with the newly forming stars 
and that their metallicities are similar. However, over time, 
the metallicity of the newly formed stellar population may 
lag behind that of the ionised gas as the gas is continually 
enriched. This implies that for younger populations with rel- 
atively large ages, the O/H abundance from Tremonti et al. 
(2004) may overestimate the true abundance of the stars. 
Even so, placing strict limits on the stellar metallicity is cru- 
cial to providing reasonable age estimates for the younger 
stellar populations and the Tremonti et al. (2004) values 
are the only independent metallicity estimates available. We 
have therefore chosen to constrain the metallicity for the 
younger stellar population to be within the 68% confidence 
interval of the O/H value given by Tremonti et al. (2004), 
and caution the reader that for galaxies with relatively large 
age estimates for their younger populations, the O /H abun- 
dances may be somewhat inflated. 

For the older population, we do not have a comparable 
constraint for the stellar metallicity. While accurately es- 
timating the age of the older population is not crucial to 
our main goal, we may partially constrain the older popula- 
tion's metallicity using the spectral index [MgFe]' (Thomas 
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ct al. 2003). While not independent of age, [MgFe]' is more 
strongly affected by metallicity and has the added advantage 
of being nearly independent of the abundance of a-elements. 
The latter is relevant since many early-type galaxies are 
known to be a-enhanced relative to the sun (e.g., Thomas et 
al. 2003) and the Bruzual & Chariot (2003) models assume 
solar values for the abundances of heavier elements relative 
to one another. According to the Bruzual & Chariot (2003) 
models, an instantaneous burst population that is older than 
~1 Gyr will have a larger value of [MgFe]' than an actively 
forming population, regardless of the new population's age. 
This implies that the measured value of [MgFe]' for the ac- 
tual galaxy spectrum provides a lower limit for [MgFe]' for 
the older population. The Bruzual & Chariot (2003) models 
also imply that [MgFe]' will not exceed 5 A. We therefore 
constrain the value of [MgFe]' for the older population to be 
< 5 A and to be greater than [MgFe]o fts — la where [MgFc]o fts 
and a are the measured values for [MgFe]' and their uncer- 
tainties taken from the SDSS spectroscopic pipeline. 

(ii) Dust extinction: Using the observed ratio of Ha to 
H/3, we can obtain a good estimate of the internal dust ex- 
tinction suffered by the ionised gas which can be used to 
constrain the amount of internal dust extinction. To esti- 
mate the extinction using the Balmer decrement, we must 
assume an intrinsic value for Ha/H/3 which is dependent 
upon temperature. To determine the best temperature to 
use, we have estimated the electron temperature using the 
[Oil], [OIII], and H/3 emission lines according to Pilyugin 
(2001) for each SFE E galaxy with > 3a detections of all 
of these emission lines. The median temperature for these 
galaxies is 5500 K, for which the ratio of Ha to H/3 is 3.02 
for case B recombination according to Osterbrock & Ferland 
(2006). We used this assumed intrinsic ratio to compute the 
colour excess for the ionised gas, E(B — V) ga s, for each SFE 
E galaxy using the extinction curve of O'Donnell (1994) 
and Rv = 3.1. We have also used the errors in the Ha and 
H/3 fluxes from Tremonti et al. (2004) to compute the la 
uncertainty in E(B — V) ga s and constrained the fitted val- 
ues for this quantity to be within ±ler of the observed value. 
When computing the model spectra, the dust extinction was 
applied to the model stellar monochromatic fluxes using the 
extinction "law" of Calzetti (2001). This extinction law uses 
the quantity E(B — V) ga s and takes into account the wave- 
length dependence of the dust extinction as well as scatter- 
ing and the geometrical distribution of the dust relative to 
the stars. 

(iii) Stellar mass and age: While we have no reliable esti- 
mates for the stellar masses of the older and younger stellar 
populations within each SFE E galaxy, we do have estimates 
for the total masses from Kauffmann et al. (2003b). How- 
ever, it should be noted that Kauffmann et al. (2003b) mea- 
sured stellar masses using linear combinations of Bruzual & 
Chariot (2003) model spectra, implying that they do not 
provide stellar mass measurements that are independent of 
our own fitting procedure. We therefore used the Kauffmann 
ct al. (2003b) value as an initial guess for the total stellar 
mass of each galaxy, but allowed the fitted value to vary 
by ±1 dex from this value. The fraction of the total stellar 
mass attributed to the younger population, f y , was used to 
quantify the relative contribution of each population to the 
total stellar mass and was only constrained to be < 1 and 
> 0. The age of the younger population was constrained to 



be older than 1 Myr while the age of the older population 
was constrained to older than 1 Gyr; both populations were 
not allowed to be older than 13 Gyr. 

4.3 Execution and reliability 

We have compiled a routine which will interpolate among 
the suite of instantaneous burst model spectra and one of the 
suites of exponentially decreasing model spectra described 
above to produce a composite spectrum for given values of 
eight parameters: (1) the age of the older population, (2) 
the age of younger population, t y , (3) the oxygen abundance 
of the younger population, O/H, (4) [MgFe]' for the older 
population, (5) the total stellar mass, M», (6) the fraction 
of the stellar mass contributed by the younger population, 
f y , (7) the dust extinction/reddening for the ionised gas, 
E(B — V) g as, and (8) the stellar velocity dispersion, a v . Since 
we keep the velocity dispersion fixed at the value measured 
for the galaxy spectrum, this leaves seven free parameters. 
We have used the IDL routine mpfit written by C. Mar- 
qwardt to determine the best fitting combination of these 
seven parameters for each SFE E galaxy that had an O/H 
measurement from Tremonti et al. (2004), 278 galaxies in 
all. This routine uses a standard Levenberg-Marquardt algo- 
rithm to minimise x 2 , but allows for constraints to be placed 
on the free parameters and numerically computes the partial 
derivatives of the function used in the fitting with respect to 
the free parameters. For each spectrum, regions of the spec- 
trum with emission lines were not used in the fit and the 
spectrum was corrected for Galactic extinction using the r- 
band extinction from Schlegel et al. (1998), the extinction 
curve of O'Donnell (1994), and R v = 3.1. For each fit, the 
free parameters were constrained as discussed above. Exam- 
ples of the SDSS spectra and the model fits for each value of 
r can be seen in Fig. 8. Overall, the models appear to fit the 
data remarkably well (within about 0.01 dex on average), 
irrespective of what e-folding time was used. 

Since the Ha fluxes were not used to constrain the spec- 
tral fits, they provide a good opportunity to test how well 
these fits represent the SFE E galaxies. Among other spec- 
tral properties, the Bruzual & Chariot (2003) models pro- 
vide the the flux of hydrogen ionising photons, Q(H°), at 
each time step for a given model, and we have used these 
values to determine Q(H°) using the fitted ages, metallic- 
ities, and masses of the older and young populations for 
each galaxy and each value of r. According to Osterbrock 
& Ferland (2006), for a temperature of 5500 K (see §4.2) 
and case B recombination, Lho = 1-46 x 10~ 12 Q(-ff°) ergs 
s _1 . We have used this equation with the fitted values of 
E(B — V) g as to compute the model predictions for the ob- 
served Ha fluxes. In Fig. 9, we have displayed the distri- 
bution of the ratios of the model-predicted fluxes to the 
observed fluxes for each value of r. Gaussian fits to these 
distributions estimate the means to be -0.28 -0.13 -0.057 
and 0.00031 dex and the standard deviations to be 0.27, 
0.19, 0.18, and 0.17 dex for r =100 Myr, 300 Myr, 1 Gyr, 
and oo, respectively. By far, the constant SFR (r = oo) 
models provide the best statistical match to the observed 
Ha fluxes with ratios of model-predicted to observed fluxes 
that range from roughly 0.7 to 1.5 with a median near unity. 
We note that this result somewhat contradicts other simi- 
lar studies which estimate that the typical e-folding time 
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Figure 8. Examples of the model spectra fits described in §4.2. For each galaxy, we have plotted the SDSS spectrum (solid black) and 
the best fitting model spectrum (solid red) for SFR e-folding times of r =100 Myr, 300 Myr, 1 Gyr, and oo (i.e., constant SFR) in each 



column from left to right. The monochromatic luminosity is plotted in units of Lq 
(Q m ,n A ,/i)=(0.3,0.7,0.7). 



: 3.826 x 10 33 ergs s _1 and was computed assuming 



for galaxies similar to our SFE E galaxies is between 50 and 
150 Myr (Kaviraj et al. 2007; Schawinski et al. 2007). How- 
ever, we also note that while this result implies that rather 
large SFR e-folding times may be common among our SFE 
E galaxies, significant portions of the distributions in Fig. 9 
for smaller values of r occupy the same region as the distri- 
bution for r = oo. In particular, the area of overlap between 
the t — oo distribution and that for each of the other values 



of r corresponds to 59%, 87%, and 96% of the total area for 
t =100 Myr, 300 Myr, and 1 Gyr, respectively. This implies 
that while our fitting results are generally consistent with 
our SFE E galaxies having a tendency for large values of r, 
we cannot rule out shorter e-folding times for many of the 
galaxies. 

To test the assumption that the limits placed on the 
fitted parameters lessened or eliminated the possible degen- 
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Figure 9. For the 278 SFE E galaxies with gas-phase O/H values 
from Tremonti et al. (2004), the distribution for the ratio of the 
the Ha flux predicted by each model fit assuming an electron 
temperature of 5500 K and case B recombination (Osterbrock & 
Ferland 2006) to the observed Ha flux. The distributions for each 
SFR c-folding time, r, arc plotted separately with Gaussian fits 
to each plotted as dotted lines. The means from these fits are 
-0.28 -0.13 -0.057 and 0.00031 dcx and the standard deviations 
are 0.27, 0.19, 0.18, and 0.17 dcx for r =100 Myr, 300 Myr, 1 
Gyr, and oo, respectively. 



eracies between the age of the younger population, t y , and 
other parameters, we have computed the Spearman rank 
correlation coefficient among all modelled galaxies for t y and 
each of the other free parameters. Of particular concern are 
the amount of dust extinction, E(B — V) ga s, and the two 
metallicity indicators, O/H for the younger population and 
[MgFe]' for the older population. This is because increases 
in any of these parameters will in general make the spec- 
tra redder, which can roughly mimic the effects of making 
the populations older. Therefore, if the fitting process has 
incorrectly inflated the values of these parameters and com- 
pensated by artificially lowering the fitted ages, we would ex- 
pect to see anti-corrlations between t y and these parameters. 
There is a small correlation between t y and E(B — V) ga s for 
the model fits with r < 1 Gyr with correlation coefficients of 
0.10, 0.26, and 0.15 for r =100, 300, and 1000 Myr, respec- 
tively. For the constant SFR model fits, there is a negligible 
anti-correlation with a coefficient of about -0.017. Both O/H 
and [MgFe]' are somewhat correlated with t y with correla- 
tion coefficients for O/H of about 0.3 for r =100 and 300 
Myr and about 0.1 for r =1 Gyr and oo and coefficients for 
[MgFe]' of about 0.3 for all star formation histories. The rel- 
atively small magnitudes of these correlations, coupled with 
the fact that they are positive, implies that the constraints 
placed on these parameters have done a reasonable job of 
eliminating problematic degeneracies between t y and metal- 
licity or the amount of dust extinction. 

A possible correlation between the age of the younger 



population, t y , and the fraction of the stellar mass attributed 
to the younger population, f y , is a special concern since since 
such a correlation can exist for both physical and artificial 
reasons. For instance, one would naively expect galaxies that 
have been forming stars longer to have larger values of both 
t y and f y , on average. On the other hand, as the newly 
formed population ages, its mass-to-light ratio will increase, 
implying that if the age of the population is over or under 
estimated, the fitting routine may be able to compensate by 
increasing or decreasing f y , introducing an artificial correla- 
tion. To explore this, we have used the covariance matrix for 
the seven free parameters computed by mpfit to compute the 
correlation coefficient for t y and f y for each spectral fit. The 
mean coefficient is about 0.7 for all star formation histories, 
indicating a significant but moderate correlation between 
these parameters which may be mostly artificial. The effect 
of this correlation on the model results is discussed below. 

4.4 Model results 

The main focus of the spectral fitting procedure detailed 
above was to estimate the distribution of ages for the 
younger stellar populations within the SFE E galaxies. Since 
the possible artificial correlation between t y and f y will sig- 
nificantly influence the interpretation of this main result, 
we have displayed the two-dimensional t y ,f y distribution 
among all SFE E galaxies for each value of r in Fig. 10. 
From these plots, it is evident that the typical age of the 
younger stellar population increases with increasing SFR e- 
folding time, t, as one might expect. For all star formation 
histories, there is a paucity of galaxies with t y > 1 Gyr. How- 
ever, inspection of the two-dimensional histograms shows 
that for larger values of r, this may be influenced somewhat 
by the t y J y correlation discussed above. In particular, those 
galaxies with log t y /l Gyr ~ — 0.3 which have f y values no- 
ticeably larger than the location of the peak of the distri- 
bution (log f y > — 1.3) may in reality have larger values of 
t y , but the fitting routine artificially increased f y instead of 
appropriately increasing t y to fit their stellar continua. We 
note, however, that only about 8-10% of the galaxies have 
-0.4 < log ty/l Gyr < -0.2 and log f y > -1.3 for r > 1 
Gyr (see Fig. 10). This implies that while there is legitimate 
reason to doubt the estimates for t y for a relatively low num- 
ber of individual galaxies, the dearth of galaxies with t y > 1 
Gyr appears to be a valid result. 

To demonstrate that our fitting routine has indeed suc- 
cessfully separated the stellar continua of our SFE E galax- 
ies into older and younger components, we have displayed 
the distributions for the ages of both the younger and older 
populations in Fig. 11 for each value of r. For all values 
of r, the two distributions are noticeably distinct with the 
difference in ages for the two populations being about 4-5 
Gyr, on average. For r > 1 Gyr, the fraction of galaxies 
with older population ages < 2 Gyr, 3-5%, is larger than 
the < 1% of galaxies with older population ages < 2 Gyr 
found for r < 300 Myr. This is likely an indirect result of 
the consequences of the t y ,f y correlation discussed above as 
the model may force the older populations in some galaxies 
to be younger to compensate for underestimating t y . How- 
ever, as noted above, this is only relevant for a relatively 
low number of SFE E galaxies and does not invalidate the 
conclusion that there are relatively few SFE E galaxies with 
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Figure 10. The two-dimensional histograms for t y and f y (see 
§4.4) for the SFE E galaxies for each SFR e-folding time, t. The 
Spearman rank correlation coefficient for t y and f y for each value 
of t is printed in each panel as well. 



Figure 11. The distributions for the ages of the younger (black 
line) and older (grey line) stellar population derived from the 
model fits to the stellar continua (see §4) for each value of the 
SFR e-folding time, r. 



younger populations older than 1 Gyr. This main result im- 
plies that the star formation within these galaxies typically 
lasts less than ~1 Gyr. 

To demonstrate that their relatively short star forma- 
tion time scales will assure that the SFE E galaxies will 
become K+A galaxies, we have constructed model spectra 
using the Bruzual & Chariot (2003) models for a 1 Gyr 
episode of star formation at a metallicity of Z=0.05 for all 
values of r. We have chosen 1 Gyr because while the typical 
younger stellar population age is ~600 Myr, the relatively 
sharp drop in the distribution of ages at ~1 Gyr (see Fig. 10 
and 11) indicates that this is a better choice for the typical 
star formation episode duration. The r = 100 Myr models 
have SFRs that decrease quickly enough that they already 
have W(Ha) < 1.6 A for ages > 1 Gyr regardless of the 
magnitude of the star formation episode (see below). This 
implies that no 1 Gyr SFR cut-off is needed for these models 
to be able to explain the lack of SFE E galaxies with t y > 1 
Gyr (see Fig. 4), and therefore, no such cut-off was imposed. 
For the other values of r, such a cut-off is required and was 
produced by forcing the model SFRs to abruptly go to zero 
at a time step of 1 Gyr. 

We have combined these model spectra with the spec- 
trum of a 5 Gyr (the median age of the older populations 
from the model fits) instantaneous burst population with 
Z=0.05 for episodes of several different strengths which we 
characterise by the value of f y at 1 Gyr (i.e., the total frac- 
tional increase of the stellar mass) . We have measured both 
W(Ha) and HSa for these composite spectra, assuming 
T e = 5500 K and case B recombination to convert the val- 
ues of Q(H°) into Ha luminosities. The paths galaxies would 
take in the HSa, W(Ha) plane according to these model val- 
ues are shown in Fig. 12 along with the measured data for 
the SFE E galaxies for each value of r. Each model track 



starts at HSa = —2.5 A and W(Ha) = 1.8 A, and loops 
clockwise around the H5a,W(Hoi) plane, with larger loops 
corresponding to larger values of f y . These model tracks 
show that only those galaxies whose stellar mass will be in- 
creased by about 2-5% or more, depending on the value of r, 
will become K+A galaxies according to our definition (i.e., 
the model tracks will pass through the shaded area in Fig. 
12). According to the fitted model spectra, the median frac- 
tional increase in the stellar mass for our SFE E galaxies for 
1 Gyr of star formation is ~ 4% with about 78% having a 
fractional increase of more than 2%. This is consistent with 
what has been found by similar, complementary studies con- 
ducted with early-type, SDSS galaxies (Kaviraj et al. 2007; 
Nolan et al. 2007; Schawinski et al. 2007). It should be 
noted that this only applies to the stellar mass within the 
SDSS fibre aperture since we do not have constraints on the 
full spatial extent of the star forming regions within these 
galaxies. The model tracks also show that the K+A phase 
ends significantly less than 1 Gyr after the end of the star for- 
mation episode, depending of the strength of the episode and 
the SFR e-folding time. This is unexpected given the fact 
that our SFE E and K+A E galaxies are found in roughly 
the same numbers. However, we note that for this exercise, 
we have assumed that the SFR instantaneously goes to zero 
for t > 300 Myr after 1 Gyr of star formation, when in real- 
ity, the SFR may decrease more gradually. This implies that 
the time a galaxy will spend as a K+A galaxy estimated 
from the panels in Fig. 12 is a lower limit. We also note that 
the strict emission line ratio (see §2.1 and Fig. 1) and mor- 
phological (sec §2.2) selection criteria used to identify SFE 
E galaxies may have caused a significant underestimation of 
the number of E and SO galaxies with recent star formation. 
This implies that this result does not necessarily contradict 
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Figure 12. For the SFE E galaxies, we plot W(Hce) versus H5a- Here, W(Ha) was corrected for extinction according to log 
W(Ha)/W(Ha) i, a = 0AE(B — V) (Calzctti 2001) where the colour excess was derived from the Balmer decrement (see §4.2). Also 
plotted are model tracks for a 1 Gyr episode of star formation added to a 5 Gyr old stellar population for different values of f y at 1 
Gyr (i.e., the total, fractional increase of the stellar mass caused by the star formation episode; see §4.3) for each SFR e-folding time, 
t. For t > 300 Myr, the SFR was forced to instantaneously go to zero after 1 Gyr (see §4.4). Times step at 1 Myr, 10 Myr, 100 Myr, 1 
Gyr, 1.5 Gyr, and 2 Gyr arc highlight with different point styles along each track. Our definition of K+A galaxies is indicated by the 
shaded region with two upper limits, the original upper limit given in Fig. 5, and the original limit with an extinction correction applied 
according to Calzetti (2001) using the median E(B — V) derived from the Balmer decrement for the SFE E galaxies (see §4.2). 



the numbers SFE E and K+A E identified for this study 
(see Table 1). 



5 STAR FORMATION TIME SCALE AND THE 
POSSIBLE FORMATION OF NUCLEAR 
STELLAR DISKS 

On average, the star formation episodes within our SFE E 
galaxies will only increase their stellar masses by roughly 
4%. The SFE E galaxies have a median stellar mass within 
the area of the SDSS 3 arcsec fibre aperture of about 6 x 10 9 
M Q , implying that once the star formation has ceased within 



these galaxies, the mass of newly formed stars will be about 
2x10 s M©. For reference, for the SFE E galaxies, the radius 
of the SDSS fibre aperture ranges from 0.1R e to 2R e , where 
R e is the half-light radius from the Blanton et al. (2003b) 
Sersic fits, or 0.2 to 3 kpc with median values of 0.4R e and 
1.5 kpc. While such a relatively small increase in mass will 
do little to change the galaxies' overall appearances, the ap- 
proximate stellar mass of newly formed stars is very similar 
to what has been found for nuclear stellar disks (Scorza & 
van den Bosch 1998; Pizzella et al. 2002). The discovery 
of nuclear stellar disks in the process of being formed would 
be extremely useful as it is not known whether such disks 
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formed with their host galaxies or whether they formed out 
of gas deposited at a later epoch. Many early-type galaxies 
have been found to contain central gas and/or dust disks 
(Kormendy & Djorgovski 1989). However, only two galax- 
ies, NGC 5845 (Kormendy et al. 1994) and NGC 4468A 
(Kormendy et al. 2005) , that have nuclear stellar and dust 
disks are known. There is some indication from the colours of 
these two disks that they are younger than the surrounding 
stellar population and that they may have formed more re- 
cently. However, if the relatively weak star formation within 
our sub-sample of SFE E galaxies can be shown to be con- 
fined to compact nuclear disks, this would provide concrete 
evidence that similar disks typically formed more recently 
than the rest of the stars within their host galaxies. 

There is other evidence that the SFE E galaxies may 
harbour nuclear stellar disks that are currently forming. 
Central stellar disks are preferentially found within lower- 
mass early- type galaxies that have steep or "power-law" 
surface brightness profiles within radii of less than 100 pc 
(Lauer et al. 1995). This is in contrast to higher-mass early- 
type galaxies that have nearly flat or "core" inner profiles 
(Faber et al. 1997). In Fig. 13, we have plotted z-band 
absolute magnitude versus velocity dispersion for all early- 
type galaxies with m r <16 and for our SFE E galaxies. For 
reference, we have also highlighted the ranges in velocity 
dispersion inhabited by power-law and core galaxies taken 
from Faber et al. (1997). From this plot, it can be seen that 
there is relatively little overlap between core and power-law 
galaxies and that the SFE E galaxies lie entirely within the 
range inhabited by power-law galaxies which tend to be the 
hosts of nuclear stellar disks. 

The gas kinematics of the SFE E galaxies are also con- 
sistent with the notion that the current star formation in 
these galaxies is occurring within disk-like structures. If the 
star formation is taking place within disks, they are likely 
rotationally supported. If this is the case, the effects of in- 
clination will contribute significantly to the shape of the 
distribution of the ratio of the observed velocity width of 
the nebular emission lines, o\,(gas), and the LOS velocity 
dispersion of the stars, a v (stars). This is because the mea- 
sured stellar velocity dispersion is likely dominated by the 
older stellar population since, according to the model re- 
sults, it typically makes up about 98% of the total stellar 
mass and about 70% of the total V-band luminosity of each 
SFE E galaxy, while the emission line gas is likely only asso- 
ciated with the younger stellar population. Because of this, 
if the younger population is forming within a rotating disk, 
galaxies that have disks with smaller inclination angles will 
skew the distribution of a v (gas) ja v (stars) toward lower val- 
ues. In Fig. 14, we plot the distribution of a v (gas)/o\, (stars) 
for the SFE E galaxies, where values for o\,(gas) were taken 
from Tremonti et al. (2004) and the values of a v (stars) were 
measured using vdispfit (see §3). From this figure, it can be 
seen that the distribution is indeed asymmetric with a larger 
fraction of galaxies having a v (gas)/a v (stars) <1. 

To explore the possibility that this asymmetry is in- 
dicative of rotating gas, we have constructed the following 
simple models. We first made the simplifying assumption 
that the stellar distribution is "pressure" supported, i.e., 
a v (stars) = v c , where v c is the circular velocity. We then 
allowed the gas to have both random and ordered (i.e., ro- 
tational) velocity components, which are related to the cir- 




Figure 13. The absolute z-band magnitude versus the stellar ve- 
locity dispersion for all early-type galaxies (black points) and for 
SFE E galaxies (white circles). The ranges in velocity dispersion 
inhabited by power-law and core early-type galaxies from Faber 
et al. (1997) are marked for reference. 



cular velocity according to vl = a% + Vq where <j and vo 
are the random and ordered components, respectively. In 
this case, the LOS velocity dispersion of the gas is given 
by a„(gas) = ^/o-Q + UQsin(i) 2 where i is the inclination an- 
gle for the axis of rotation of the gas. Combining the above 
relations yields the following 



Mggs) V 
a v (stars) J 



fesin(i)] 2 + l 



(2) 



Using equation (2), we have constructed several model 
distributions for a v (gas)/ a v (stars) for different values for 
Vol g . For all models, we assumed that i was distributed 
according to sin(i), the distribution one would expect for 
random orientations (Hubble 1926). To incorporate the ef- 
fects of observational uncertainty, we fit a Gaussian to the 
distribution of errors in a v (gas)/<r„ (stars) for the SFE E 
galaxies and convolved each model distribution with this 
Gaussian. These model distributions are plotted with the 
data in Fig. 14 with the Gaussian kernel plotted as a model 
where v /<j = 0. The figure demonstrates that in order to 
reproduce the observed distribution, some amount of rota- 
tion is required for the gas. The model curves show that a 
single value of v /a cannot reproduce the entire observed 
distribution and that values ranging from about two to ten 
are required. To better illustrate this, we have also included 
another model curve (the red curve in Fig. 14) which used 
a log- normal distribution of values for v /a with a me- 
dian of two and a dispersion of 0.2 dex. This model curve 
is a close approximation to the observed distribution for 
0.2 < a v (gas)/<7„ (stars) < 1.3 and implies that v /a ranges 
from 0.8 to five with a typical value near two. 

We also constructed model distributions for a scenario 
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Figure 14. The distribution of the ratio of the LOS velocity 
dispersion of the ionised gas to that of the stars for the SFE E 
galaxies. The curves represent model predictions computed us- 
ing equation (2) for different values of vq/oq for the gas, where 
vo and <To are the rotational and random velocity components, 
respectively (see §5). The red curve represents a model that as- 
sumes a log- normal distribution for vo/ao with a median of two 
and a standard deviation of 0.2 dex. 



where the stars also have a rotational component and v /a 
for the gas has a log-normal distribution. For this exercise, 
we fixed v /<j for the stars at 0.65, the amount of rota- 
tion needed to produce the maximum allowed ellipticity of 
0.3 (see §2.2). To reproduce the red curve in Fig. 14 with 
this model, it was necessary to increase the median value 
of Vo/uo for the gas to 2.3-2.5, depending on how the ro- 
tation axes of the stars and gas were aligned. This implies 
that the estimated typical value of v /a ~ 2 for the gas 
is a lower limit. This further demonstrates that not only 
does the observed distribution indicate the presence of ro- 
tation within the gas, it also implies that the rotational 
component of the gas velocity is the dominant component. 
We note that the observed distribution does have a tail for 
<j v (gas) /<7„ (stars) > 1.6 that is not reproduced by any of the 
models. However, this tail consists of only 2% (7 galaxies) 
of the SFE E sub-sample, all but one of which have errors 
in their gas velocities that are unusually large (between 15% 
and 100%). We therefore conclude that the data are consis- 
tent with the star forming regions in SFE E galaxies typi- 
cally being contained within rotationally supported disks. 

Since our SFE E galaxies have no prominent disks vis- 
ible in their g-band images, it is reasonable to infer that 
the disks inferred from Fig. 14 may be relatively small and 
similar to known nuclear stellar and gas disks. In fact, the 
estimated typical value of v /<j ~ 2 for the gas is consistent 
with what has been observed for the stellar populations in 
the centres (radii < 100 pc) of early-type galaxies with iden- 
tified nuclear stellar disks where v /a ~ 1 (see, e.g., Kra- 
jnovic & Jaffe 2004). However, the fact that the observed 



distribution of cr v (gas)/er„ (stars) peaks near unity implies 
that the ionised gas probes the same range in circular ve- 
locity as the stars. Since the SDSS fibre aperture extends to 
radii of at least 200 pc and up to 3 kpc (or 0.1 to two effec- 
tive radii) for our SFE E galaxies, that result seems to imply 
that the disks may be somewhat extended. But, we must also 
recognise the fact that <r„(gas) and a v (stars) are both lumi- 
nosity weighted velocity dispersions. Since the mean Sersic 
index for our SFE E galaxies is about 4.5, the stellar velocity 
dispersions must be biassed toward smaller radii, implying 
that the Ha surface brightness profiles must be similarly bi- 
assed. Therefore, even though some spatially extended star 
formation may exist within these galaxies, they likely have 
some centrally concentrated star formation. This makes our 
SFE E galaxies excellent candidates for follow-up narrow- 
band imaging and spectroscopy to search for nuclear disks 
with active star formation. 



6 DISCUSSION AND CONCLUSIONS 

Using the vast amount of galaxy imaging and spectroscopic 
data available from the SDSS, we have selected a relatively 
large, clean sample of actively star-forming E and SO galax- 
ies. These galaxies make up about 3% of all early-type SDSS 
galaxies with m r <16 and 0.7% of all Main (Strauss et al. 
2002) galaxies with m r <16. These galaxies appear to be less 
massive on average than typical early-type galaxies which is 
consistent with what has been found by the similar, com- 
plementary study of Schawinski et al. (2007). This is also 
consistent with what has been found by other authors (e.g., 
Treu et al. 2005; Graves et al. 2007), that the fraction of 
stellar mass formed at low redshift is substantially higher for 
lower mass early-type galaxies. The number of star forming 
E and SO galaxies is also comparable with the number of 
morphologically similar K+A galaxies. The relatively com- 
plex selection criteria used to select both star forming E and 
SO galaxies and K+A galaxies (see §2) make it difficult to 
compare the sample sizes in a statistically robust way. How- 
ever, the fact that the numbers of both types of galaxies are 
not vastly different does allow for the possibility that the 
star forming E and SO galaxies are the progenitors of E and 
SO K+A galaxies. 

Through modelling of the stellar continua of our SFE 
E galaxies, we have demonstrated that this is indeed plausi- 
ble. We have shown that our star-forming E and SO galaxies 
likely form stars for ^ 1 Gyr. The model results demonstrate 
that the factor that most strongly determines whether or 
not one of these galaxies will become a K+A galaxy is the 
magnitude of its star formation episode. According to our 
K+A definition, only star formation episodes that increase 
the stellar mass by roughly 2-5% or more will lead to K+A 
galaxies, which appears to be the case for about 80% of the 
SFE E galaxies. This suggests that in terms of star forma- 
tion episode duration, all star-forming E and SO galaxies 
are similar to E and SO K+A galaxies, even though some of 
them will not have Balmer absorption lines that are strong 
enough for them to be classified as K+A galaxies after their 
star formation has ceased. 

With the implication that our samples of star-forming 
and K+A E and SO galaxies are linked in an evolutionary se- 
quence comes a large potential for followup observations. For 
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instance, as discussed in §5, there is evidence to suggest that 
star-forming E and SO galaxies may be the sites of nuclear 
stellar disk formation. High spatial resolution narrow-band 
images of the star-forming E and SO galaxies will be able to 
determine whether or not this is true. If so, E and SO K+A 
galaxies will be excellent candidate hosts of nuclear stellar 
disks. If such disks are responsible for the galaxies' strong 
Balmer absorption, then they should appear relatively young 
(i.e., bluer than the surrounding stars); this may also be ex- 
plored will high spatial resolution followup imaging. 

Identifying relatively large samples of SFE E galaxies 
with these types of disks will prove extremely useful for mea- 
suring central black hole masses with future spectroscopic 
observations (e.g., using NIRSpec with the James Webb 
Space Telescope) for relatively low mass E and SO galaxies. 
This is because the kinematics of these nuclear disk struc- 
tures are dominated by rotation and velocity anisotropy is 
relatively unimportant (Kormendy et al. 2005). However, 
estimates for BH masses using nuclear stellar disks are based 
on luminosity weighted stellar velocities measured from op- 
tical spectra which are affected by radial variations in the 
stellar population (Kormendy & Richstone 1995). Rotation 
velocities measured for gas using emission lines are not af- 
fected by this phenomenon. In this sense, the interpretation 
of the kinematics of any emission line gas associated with a 
nuclear disk in regards to the mass of the central BH will 
be much more straightforward. Identifying SFE E galaxies 
that contain such disks will therefore provide a sample for 
which future spectroscopy can be used to reliably estimate 
central BH masses. 
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